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The interplay between antiferromagnetism and d-wave superconductivity is studied in a mean- 
field approximation for a generic microscopic Hamiltonian with short-range repulsion and near- 
neighbor attraction. In the presence of competing microscopic interactions, the phase boundaries 
of antiferromagnetic and superconducting states are significantly modified in some region of the 
doping-temperature plane. The transition between superconductivity and antiferromagnetism oc- 
curs through a phase where both order parameters coexist with a third, dynamically generated, 
spin-triplet amplitude. This dynamical generation of a new order parameter is not restricted to a 
system with antiferromagnetism and d-wave superconductivity, but is a generic feature for fermionic 
systems. The dynamically generated spin-triplet order parameter is found to be robust to variations 
in the mean- field Hamiltonian. 
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The close interplay between antiferromagnetism (AF) 
and superconductivity (SC) has been of great interest 
both to experimentalists and theoretical physicists for 
many years. The competition of two different order pa- 
rameters has been observed in various systems such as or- 
ganic superconductors, |^ ternary compounds, Q heavy 
fermion compounds ^ as well as high temperature su- 
perconductors 1^ in the latest. For some compounds, 
such as the cuprates, the AF and SC phases are close 
to each other but separated, while for other compounds 
such as organic superconductors, ternary compounds and 
heavy fermion compounds the two phases often touch and 
sometimes even coexist in the phase diagram. This com- 
petition problem of two different orders was extensively 
studied in a mean-field approximation first for s-wave 
superconductors Q and later for anisotropic supercon- 
ductors 1^,0. Recently the same problem including the 
fluctuation effect of AF and SC was considered by On- 
oda and Imada [^j9| in the framework of self-consistent 
renormalization theory. 

A few years ago Zhang and his coworkers pro- 
posed that the SC and AF phases, particularly, of the 
high Tc cuprates might be related to each other in a uni- 
fied theory based on S0(5) symmetry. In this S0(5) 
theory, the SC components of the five dimensional super- 
spin vector are transformed to the AF components by, so 
called, 77 operators. These operators carry charge 2 and 
spin 1, and total momentum Q — (tt, tt) in two dimen- 
sions. Further the above authors showed that a sharp res- 
onance of 41 meV observed in inelastic neutron scattering 
experiments [ |T2| might be a direct consequence of the ex- 
istence of TT operators. To the author's knowledge, the 
possibility of spin-triplet amplitude with total momen- 
tum Q often called 7r-triplet in the literature, was first 



found by Psaltakis and Teuton in their mean-field 
study of the competition problem of AF and SC. These 
authors showed that when a spin-density wave (SDW) 
coexists with spin-singlet Cooper-pair superconductivity 
among the same electrons, a non-zero spin-triplet am- 
plitude must appear. Recently a similar problem was 
studied by Murakami and Fukuyama by invoking 

the g-ology diagrams in one-dimensional electron system. 
But in the above two works the authors did not study 
the condition of appearance of the spin-triplet amplitude 
and its robustness. Furthermore in previous mean- field 
studies the manner how the AF-only and SC-only phase 
boundaries are modified into new phase boundaries in 
the presence of two competing interactions, was not well 
investigated. Thus the purpose of this paper is twofold: 
First, detailed phase diagrams are presented in the dop- 
ing (x = 1 — n)-temperature (T) plane in the presence 
of competing interactions and second the condition of 
appearance of the spin-triplet order parameter and its 
robustness are clarified. 

In this paper we study a phenomenological mean-field 
Hamiltonian where antiferromagnetism and d-wave su- 
perconductivity are built in on equal footing. 

H - E(^£ - ^'K.%^ + c^E«T^^T«iCa) + H.c.) 

- FE(At,(A<i,.) + H.C.) - I^E(^m(^m) + H.c; 



(1) 



A possible spin-triplet pair amplitude {^t,i) is also con- 
sidered in general for self-consistency and its significance 
will become clear later. In particular, we will see that it 
is dynamically generated, even when W = Q. The sin- 
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glet Ad^i and triplet At^^ pair destruction operators are 
defined as 



where 



(2) 



and the structure factor g{S) is chosen to have a d-wave- 
hke form such as 



1/2 if J =(±1,0), 
1(6) = { -1/2 ifj=(0,±l), 
if otherwise. 



(3) 



Later it will be also clear why the same structure factor 
g{S) is introduced in the spin-triplet term. The strength 
of AF, SC, and spin-triplet pair interactions is governed 
by U, V, and W, respectively, which are all positive 
in our study. Essentially the same form of the mean- 
field Hamiltonian was previously studied by Psaltakis 
and Fenton, ||l3| and Kato and Machida, |^ and Mu- 
rakami and Fukuyama ||l^ except for some slight differ- 
ences, ^j: is a tight binding energy dispersion given as 

= — 2t(cos fcj; -I- cos ky) — At' cos cos ky, where t and 
t' are hopping constants for nearest neighbors and next 
nearest neighbors. And /i is the chemical potential which 
controls the particle density n. 

In a mean-field approximation, three different order 
parameters m, s, and t corresponding to the three differ- 
ent interactions may be defined in the following way: 

{ctcr^i.<y) = ("i,"-) = + crmcos{Q ■ Ri) 
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g[5){ci+s,]Ci^l + Ci+s,iCi,]) = tcos{Q ■ Ri) , (4) 



where Q is the (commensurate) antiferromagnetic wave 
vector equal to (tt, tt) in two dimensions. Now the mean- 
field Hamiltonian Hmf is quadratic in the original elec- 
tron operator and in terms of a new four component field 
operator tp^ it becomes bilinear 



Hmf — 



k 



where 



+ _ 



The matrix AIj^ is given as 

/ Ej^ Vs(j){k) 

Vs<p{k) -e^ 
-Urn --Wt(t>{k) 
\ Wt(j){k) -Um 



(5) 



(6) 



-Um Wt(j){k) \ 

Wt(t){k) ~Um 

£fe+Q^ -Vs^ik) 

Vs(l){k) -ej:^0 ) 



^k=^k-^ 



(8) 



and 4>{k) = cos k^ — cos ky is the Fourier transform of 
g{&\ The prime symbol on the summation requires the 
summation of wave vectors in halj of the first Brillouin 
zone, in order to take into account the doubling of a 
magnetic unit cell in the presence of (commensurate) AF 
order. The constant energy shift Eq depending on to, s, 
t, and is given as 



(9) 



where N is the total number of lattice sites. The energy 
eigenvalues of A-Ij: yield four energy dispersions ±E±{k), 

E±{k) = [ (e| + £|^^)/2 + {Um)' + [Vsmf 

+ {Wtcf,{k))'±g{k) ]^^' , (10) 

where g{k) is given as 

9ik) = [ (el - + ^'k - e~,^^)'iWtmf 

+ ii^k + H+Q)iUm) + 2{Vsm){Wtm)? Y'^ ■ 

(11) 

When s = t = Oorm = i = 0, the energy eigenvalue 
reduces to that of SDW or BCS approximation of the 
corresponding interactions. 

The free energy is easily obtained either from the trace 
formula or from the Feynman theorem 

F = -2r^ ^ l°g(2cosh%^) + E, . (12) 



Now the three mean-field equations are obtained by the 
stationary condition of F with respect to the correspond- 
ing order parameters ~ ^ ~ ^ ~ ^ ^^"^ '^^^ 
more unknown constant fj, is determined by the thermo- 
dynamic relation n = The resulting four equations 
are 



(^fc ^k+Q' 

2g{k) 



ii^k + ^k+Q^iUm) 



+2{Vscb{k)){Wtcf,{k))]}^{k)tanhi^^^^ 



9(.k) 



(13) 

i^k + ^k+q) 



(7) (t/m) -t-2(Fs0(fc))(m0(fc))]|-^(fc)tanh( 

J En 



pE^ik), 



(14) 
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{Wmk)){e^-e^^ 



dependent of t, 



'igih) 



9{k) 

+2{Vs(l){k)){Wt(l){k))] I (fc) tanh( 



(15) 



a=± 



2g{k) 

+2{Vs4,{k)){Wt4,{k))] \ (fc) tanh( 



(16) 



Although essentially the same form of the mean-field 
Hamiltonian was studied by Psaltakis and Fenton, fl^ ] 
and Murakami and Fukuyama fl^ , the final mean- field 
equations are somewhat different. In Ref. ||l3| the Fermi 
surface is divided into two regions. In region one the nest- 
ing condition = —Sj:^^ is satisfied and in the other 
region it is not satisfied. Thus the mean-field equations 
are different in different regions. Even after taking this 
into account there are still some discrepancies. In Ref. 

apparently there are several erroneous terms in the 
mean-field equations. As a double check for the correct- 
ness of the above mean-field equations, we also computed 
the transformation matrix relating the original c opera- 
tor to the eigen-operator 7 and with it obtained the same 
mean-field equations. All the components of the matrix 
are given in the Appendix. 

Before plunging into numerical calculations, we point 
out several distinctive features inferred from the struc- 
ture of the mean-field equations itself. First, a given 
order parameter is affected by the other ones in a more 
direct manner than for the case where only two order pa- 
rameters, AF and SC, [D are considered. For the latter 
case the influence from the other order parameters comes 
indirectly through Ea{k) and g{k). Second, in spite of 
non-vanishing driving potentials for the three different 
order parameters, only a state with a single order param- 
eter or with three order parameters together is possible. 
For the former case the proof is obvious. Suppose two 
of the three order parameters, for example, m and s, are 
non-vanishing and then examine whether the third one 
(t) can vanish or not. Since the equations have the same 
structure with respect to the three order parameters, this 
particular choice suffices for general purpose. The right 
hand of the t equation, Eq. 15, has a term which is in- 



{Vsm){£k+euQ)iUm) 



(17) 



Thus unless this term is exactly canceled by the combi- 
nation of the other terms, the third order parameter t 
does not vanish. But perfect cancellation is impossible 
because the term Eq. 17 cannot be written as a linear 
combination of the other terms in Eq. 15. Third, despite 
the vanishing driving potential for one of the order pa- 
rameters, for example, W — 0, the corresponding order 
parameter t (in this case) can be dynamically generated. 
Since the terms containing t enter with the combination 
of Wt(p{k) in the right hand of Eq. 15, only Eq.^ sur- 
vives for W — 0. As long as both AF and SC order 
parameters have non-vanishing amplitude, namely, they 
coexist, t also has non-vanishing amplitude even with 
W = 0. 

This last feature was first discussed by Psaltakis and 
Fenton They argued that when AF and SC order 

parameters (ci^Cj:_|_^ -j-), (cg|C_^j^) have non-vanishing 
value at the same time, two electrons cr , a ^, c r 1 are 
coupled to the same electron ci through the Bogoli- 
ubov transformation, thus they are not generally inde- 
pendent but instead may form another order parameter 
{^k+Q I'^-k i)' ^Q'll^ shows in which circumstances this 
may happen. The d-wave factor (f>{k) in front of the curly 
bracket in Eq. 15 comes from the spin-triplet term and 



(/)(fc) in Eq, 17 from the spin-singlet term, thus the same 
symmetry factors, or at least those which are not orthog- 
onal to each other, are required to obtain a non- vanishing 
third order parameter. At half-filling for a nearest neigh- 
bor model (t' — 0) where e^r + £i^_^_Q = 0, however, the 
dynamical generation of another order parameter does 
not happen. Since the equations again have the same 
structure with respect to the three order parameters, m 
(s) can be dynamically generated with U = {V = 0) a,s 
well. In this paper we argue that this dynamical genera- 
tion of third order parameter in the coexistence state of 
two orders is not restricted only to the case with AF and 
d-wave SC orders, but is a generic feature for fcrmionic 
systems. With a proper choice of the symmetry factor, 
a third order parameter can be dynamically generated 
without the corresponding driving potential. For exam- 
ple, for the case of SDW and p-wave SC, the third order 
parameter is 0(fc)(cj -^c^^^ where (t>{k) is a p-wave 
form factor such as sin k^ or sin ky . This may be called 
TT-singlet in parallel to the 7r-triplet for the case of SDW 
and d-wave SC. In S0(5) theory, the 7r-triplet amplitude 
is not an order parameter. It is the generator of infinites- 
imal rotations between AF and SC order parameters. In 
the mixed AF+SC+7r-triplet phase that is discussed here, 
the TT-triplet itself acquires a non-zero expectation value. 

First we start with presenting phase diagrams of the 
nearest neighbor model {t' = 0). In this paper, t, lattice 
spacing, ?i, and ks are set to be unity, and U is fixed at 
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half of the bandwidth, 4. With varying V, the shape of 
the phase diagram changes significantly. In the present 
study we restrict ourselves to a particular case where at 
half filling the ground state is an antiferromagnetic insu- 
lator, which might be more relevant especially for high 
Tc cuprates. This requirement dictates that V is less 
than U — A and thus we consider two different values 
of V, 3 and 1.5. In Fig.^ and in the following figures, 
the dotted and dashed curves denote the phase bound- 
ary of AF and SC without the driving interaction term 
for the other order parameter. For instance, the dotted 
curve in Fig.|l|-(a) is the phase boundary of AF state with 

V = and the dashed one of SC state with U — 0. Fig.|l| 
clearly shows how the AF-only (dotted curve) and SC- 
only (dashed curve) phase boundaries change into new 
ones (solid curves) in the presence of two competing in- 
teractions. For F = 1.5 the major change of the phase 
boundary appears in the low doping region of the SC 
state. Near half-filling, due to strong local repulsion, 
the SDW state is formed throughout the lattice and it 
strongly breaks time-reversal symmetry, leading to the 
destruction of the SC state. But the AF-only phase re- 
mains virtually unchanged at this strength. In the low 
doping region a new phase appears where the AF and 
SC and 7r-triplet order parameters are non- vanishing, as 
described earlier. For V = 3.0 where the strengths of AF 
and SC interactions become comparable, both the AF- 
only and SC-only phase boundaries significantly change. 
As in y = 1.5 the SC phase boundary is most strongly 
modified near half-filling, while the AF phase boundary 
far away from half-filling where AF order is most vulnera- 
ble to other interactions. Because of the larger value of V, 
the SC phase region becomes wider than that oiV — 1.5. 
At a particular point where four different phases (para- 
magnetic metal, SC, AF, and AF-|-SC-|-7r-triplet phases) 
meet, a tetracritical point ||l6| , |ll|l appears. That tetra- 
critical point does not have SO (5) symmetry in general, 
but it has SO(3)(8)U(l) symmetry. S0(5) symmetry is 
realized only when t' = 0. jl^ 

In Fig.|^ we consider phase diagrams with the next 
nearest neighbor hopping {t' = —0.3). This value 
is close to what is suggested for YBa2Cu307_a and 
Bi2Sr2CaCu208+5. For small V the coexistence region 
is more suppressed near half-filling than that for t' = 
in Fig|^(a) because the effective density of states for SC 
is reduced. For ^ = 3, however, the coexistence region 
is more strongly modified far away from half-filling than 
in the t' = case, since the increased effective strength 
for SC in this doping range pushes the AF phase to the 
left. The slight increase of the AF region near a; = 0.3 
with t' = —0.3 is presumably due to an artifact of the 
commensurate SDW approximation. 

As in the previous figures, the dotted and dashed 
curves in Fig.^ denote the order parameters of AF with 

V = and those of SC with U — 0, respectively. The 
order parameters are plotted in Fig.0-(a) as a function of 



doping {x) with temperature fixed at 0.3. As can be in- 
ferred from Fig|l]-(b), the order parameters of AF and SC 
states are most strongly affected far away from half-filling 
and near half-filling, respectively, by the presence of com- 
peting interactions. Where two order parameters (AF, 
SC) coexist, the third 7r-triplet order parameter appears 
without its driving potential (W = 0). Its amplitude 
is maximum when those of AF and SC become compa- 
rable as is implied in Eq.^ and the maximum value is 
approximately one fifth of those of AF and SC at the cor- 
responding doping level. When one of the order param- 
eters appears or disappears, the other order parameter 
changes in its slope, a characteristic feature of a mean- 
field approximation. In Fig.||-(b) the order parameters 
are plotted as a function of temperature at fixed doping 
0.1. Because of their competing nature, both order pa- 
rameters are reduced at low temperatures from those of 
AF-only and SC-only phases. 

In Fig the phase diagram is shown for W = 3, V = 3, 
and t' = 0. As in the previous figures, the dotted, dashed, 
and long-dashed curves are phase boundaries for AF with 

V = W = 0, foT SC with U = W = and for 7r-triplet 
with U = V = 0. Due to the W term, AF and SC 
regions are increased compared with those for W = case 
(Fig.|l|-(b)). This is because the W term increases the 
stability of AF and SC by converting SC (AF) electrons 
into AF (SC) electrons instead of destroying them. One 
of the interesting features is that a 7r-triplet phase does 
not appear as a separate phase for W — 3, although 
W is as large as V. In order to create such a phase, 
W > 4.5 is required. Then a SC-only phase disappears 
in the phase diagram and the roles of SC and 7r-triplet are 
interchanged. In this situation the SC order parameter 
appears only where AF and 7r-triplet phases coexist. 

In the last figure the order parameters are plotted 
as a function of doping and temperature for W = 3, 

V — 3, t' — 0. As before the dotted, dashed, and 
long-dashed curves are the order parameters for AF with 
F = = 0, for SC with U ^0 and for 7r-triplet 
with U — V — 0, respectively. Because of the W term, 
the critical concentration of AF order at T = 0.3 (Fig.||- 
(a)) is increased from x = 0.2 for = to x = 0.27, 
much closer to that of AF order with V — W = (dotted 
curve). Thus the order parameter of AF with = 3, is 
similar to that of AF with V = W = 0. The maximum 
value of TT-triplet order parameter is 0.051, 0.064, and 
0.091 for W = 0, 1.5, and 3.0, respectively. This indi- 
cates the robustness of dynamically generated {W = 0) 
part of the 7r-triplet order parameter. 

Before closing several comments are in order. In or- 
der to find the nature of phase transition at the phase 
boundary of AF (or SC) ^ AF-|-SC-|-7r-triplet, we cal- 
culated the free energies of the three different phases 
(AF, SC, and AF+SC+7r-triplet) at constant tempera- 
ture and chemical potential. Inside the phase boundary 
the free energy of AF-|-SC-|-7r-triplet is lowest [|l8| and 
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its slope does not change as the temperature crosses the 
phase boundary. Thus in a mean-field approximation 
the phase transition of AF (or SC) AF+SC+Tr-triplet 
is second order. This is consistent with the earlier 
result of Kato and Machida. Q These authors classified 
anisotropic pairing states enumerated group theoretically 
into two categories: less-competing and more-competing 
states according to the combined symmetry of the SC 
order parameter (parity and translational symmetry de- 
termined by the nesting vector). Less-competing states 
yield second order phase transition between AF (or SC) 
and ^ AF-|-SC+7r-triplet, while more-competing states 
yield first order phase transition. D-wave SC belongs to 
the less-competing states according to their classification. 
In this paper only a commensurate magnetic structure is 
considered in the SDW approximation. From more so- 
phisticate considerations, away from half- filling and 
at low temperatures the commensurate magnetic order 
should be replaced by an incommensurate one. Lastly, 
mean-field theories neglect fluctuations completely, thus 
they are inadequate in low dimensions because there is no 
long-range order at finite temperatures in one and two di- 
mensions. This pathological feature may be avoided 
by considering a small z-axis component. Since strong 
fluctuations still remain due to quasi two dimensional 
nature, the phase boundaries will be shrunk drastically. 
For AF state the boundary will be squeezed toward the 
low temperature and low doping, while for SC state to- 
ward the low temperature. Accordingly the coexistence 
region may remain in a very small region of the phase dia- 



gram or even disappear. |22| Then the mean-field critical 
temperatures will be replaced by a crossover tempera- 
ture T* below which a pseudogap will appear in the 
single particle spectral function. This precursor effect of 
the ground state can be caused by antiferromagnetic or 
superconducting fluctuations depending on the location 
in the phase diagram. Likewise, a strongly fluctuating 
TT-triplct amplitude will manifest itself in the region sur- 
rounded by the AF and SC mean-field phase boundaries. 

In conclusion, we have studied the close interplay be- 
tween antiferromagnetism and d-wave superconductivity 
in a mean-fleld approximation. In the presence of com- 
peting interactions the phase boundaries of AF and SC 
states are significantly modified in some region of the 
doping-temperature plane. In the coexistence region of 
AF and d-wave SC a new order of spin-triplet pair is dy- 
namically generated without its driving interaction. We 
argue that this feature is not restricted only to the cur- 
rent problem, but is a generic feature for fcrmionic sys- 
tems. The dynamical generation of spin-triplet order pa- 
rameter is found to be robust to variations in the mean- 
field Hamiltonian. Fluctuation effects to the mean-field 
result are briefiy discussed. 
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APPENDIX: TRANSFORMATION MATRIX 

The transformation matrix Aij relating the original c 
operator to the eigen-operator 7 is defined by 



(Al) 



The components of the first column vector are given as 

{(E+ik) + ej:)iE+ik) + ej:^^){E+{k) - sj:^^) 
-{E+{k) - ££+^)(C/m)2 - {E+(k) + ej'){Vsmf 
-{E+ik) + ej:^^){Wtmf + 2{Um){VSm) 



A2.i= {{VSm)[iE+{k) + ej:^^){E+{k) - sj:^^) 
-{Umf - {VS(t>{k)f + {Wt(t){k)Y] 

+2{Um)iWtm)eUQ)/^T. 4i 

A3,i= -{{Um)[(E+{k) + ej:)(E+{k) + e~^^^) 
~{Umf - iVS(f>{k)f - {Wt(f>{k)f] 
+ {Vsm){Wtm){2E+{k)+e^^^ + ej:)) 

^4,1- {{Wtct>{kME+{k) + ef:){E+{k) - e,,^^) 
-{Umf + {VS(t>{k)f - {Wt(t){k)Y] 

+ [Um){Vsm){e^,^Q + e£))/^^4i , (A2) 

where the square of normalization of the vector is 
Y,Al^^2E+{k)g{k)[{E+{k)+ej:) 

i 

{2g{k) + el - el^.) + 2{Umf{sj: + e^,_^^) 

+4.iUm)iVS^ik))iWtcj){k)) + 2(Wt^{k)fie^ ~ ej:^^] 

(A3) 

By symmetry the second column vector Ai_2 is obtained 
from the first one as follows: Ai^2 — —^2,17^2,2 — 
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^i,ii^3,2 = ^4,1,^4,2 = —^3,1- The third column vec- 
tor may be expressed in terms of the first one: A1.3 = 

^3,1,^2,3 = -^4,li^3,3 = ^1,1' ^4,3 = -^2,1 with 

the following replacement in A^^i, <-> s^^^^, g{k) <-> 

—g{k), E+{k) ^ E-{k). Finally the fourth column 
vector Ai^4 is found from the third one by symmetry 

yll^4 = ^2, 37^2, 4 = —^1, 3,^3, 4 — —^4, 31^4, 4 = ^3,3- 
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FIG. 1. Phase diagrams in doping (x = 1 — n) and temper- 
ature (T) plane for (a) V = 1.5, and (b) V = 3, with C/ = 4, 
W = 0, and t' = 0. The dotted and dashed curves denote the 
phase boundaries of AF with V = and of SC with U = 0, 
respectively. 
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FIG. 3. Order parameters at (a) T — 0.3 and (b) x = 0.1, 
tor U = 4, V = 3, W = 0, and t' = 0. The dotted and dashed 
curves denote the order parameters of AF with V = and of 
SC with {7 = 0, respectively. 
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FIG. 2. Phase diagrams in doping (x = 1 — n) and temper- 
ature (T) plane for (a) V = 1.5, and (b) V = 3, with U = 4, 
W — 0, and t' = —0.3. The dotted and dashed curves de- 
note the phase boundaries of AF with V = and of SC with 
U = 0, respectively. 



FIG. 4. Phase diagram in doping {x = 1 — n) and tem- 
perature (T) plane for f/ = 4, V = 3, = 3, and t' = 0. 
The dotted, dashed, and long-dashed curves denote the phase 
boundaries of AF with V = W = 0, oi SC with U = W = 0, 
of TT-triplet with U = V = 0, respectively. 
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FIG. 5. Order parameters at (a) T = 0.3 and (b) x = 0.1, 

for [/ = 4, V = 3, W = 3, and t' = 0. The dotted, dashed, 
and long-dashed curves denote the order parameters of AF 
with y = = 0, of SC with U = W = 0, oi 7r-triplet with 
U = V = 0, respectively. 
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